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Abstract 

We consider products of two 2- manifolds such as S 2 xS 2 , embedded 
in Euclidean space and show that the corresponding 4- volume preserv- 
ing diffeomorphism algebra can be approximated by a tensor product 
SU(N) (g) SU(N) i.e. functions on a manifold are approximated by 
the Kronecker product of two SU(N) matrices. 

A regularization of the 4-sphere is also performed by constructing 
_/V 2 x N 2 matrix representations of the 4-algebra (and as a byproduct 
of the 3-algebra which makes the regularization of S 3 also possible). 

1 Motivation 

Matrix regularization of surfaces is a remarkable statement which says 
that there exists a map from functions on a surface (closed, orientable) 
to SU(N) matrices such that structure constants of the area preserv- 
ing diffeomorphisms (APD) on that surface can be approximated by 
the structure constants of SU (N) . The statement was put froward in 
early 80s by Goldstone and Hoppe [lj in the case of a sphere. Later on 
it was generalized to the case of a torus [21 [3j H] , higher genus Riemann 
surfaces [U [6] and arbitrary Kahler manifolds [8] . A very important 
property of matrix regularization is that among closed, connected and 
orientable surfaces the underlying matrix group is SU(N) regardless 
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of the topology. On the other hand, the topological information is 
encoded in the way one proceeds with the N — > oo limit therefore 
the global information about the surface is in the limit. At finite but 
large ./V one can extract topological information about the surface by 
analyzing the spectral properties of the matrices as observed in [9j. 

Originally, matrix regularization was invented in order to consis- 
tently quantize a membrane pQ . However it is clear that the framework 
itself, a priori has nothing to do with quantization: it is a mathemat- 
ical framework in which one can consider functions on surfaces as a 
limiting case of certain NxN matrices. As such it can be applied both 
to classical physics by considering the regularized classical equations 
of motion [lOl QZEJ 02] and quantum physics by regularizing the quan- 
tum theory of fields on the surface - as a result obtaining a tractable 
quantum mechanical system [H [13]. In the last few years it has been 
discovered that many aspects of the differential geometry of (embed- 
ded) Riemannian 2-manifolds can be formulated by that procedure 

nana Eg. 

Matrix description of a quantum membrane is clearly interesting 
from the particle-physics point of view. If elementary particles are 
extended objects then it is natural to consider membranes and their 
excitations as good candidates - the idea put forward by Dirac in 1962 
[T71 [TBI [19] where he introduced what is now called the Nambu-Goto 
type action for membranes and considered the possibility that leptons 
may be understood as excitations of the membrane ground state (the 
electron) . 

On the other hand, it is quite natural to search for generalizations 
of matrix-like regularization to higher-dimensional manifolds. This 
problem, although very well known [20J, is still unsolved. Moreover 
from the point of view of quantum gravity one is more interested 
in a possible regularized description of space-time itself which favors 
4-manifolds among others. If, for simplicity, one assumes that the 
4-manifold under consideration is compact and that it has Euclidean 
signature then one arrives at as simple question as: what is the reg- 
ularized description of S 2 x S 2 ? In this paper we shall prove that 
the group structure underlying this manifold is the tensor product 
SU(N) ® SU(N) i.e. a group of N 2 x N 2 matrices obtained by a 
Kronecker product of two SU(N) matrices. The result can be gen- 
eralized to arbitrary product N x A4 of two 2-manifolds N and A4 
whose regularizations are given by SU(N) matrices. 

A way to regularize the 4-sphere is also presented by constructing 
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_/V 2 x N 2 matrix representations of the coordinates of S 4 embedded 
in M 5 . We than show that the Nambu 4-bracket evaluated on the 
coordinates can be replaced by a 4-commutator of those matrices. 
This construction can then be used to prove analogous statements for 
S 3 . 

2 Matrix regularization of surfaces 

Let us start with reviewing matrix regularization for surfaces. There 
are several excellent reviews on that subject e.g. [2T1[22] - our objective 
here is to fix the notation and the conventions. 

For the case of a single sphere it is convenient to expand functions 
on S 2 in terms of spherical harmonics Y^ipi, Lp2) 

oo I 

F(<px,ip 2 ) = Y^ a L Y m( < Pl^) 
Z=0 m=-l 

where a l m are real. The corresponding structure constants cabc of the 
APD of S 2 will be given by 

{Y A ,Y B } = c ABC Yc 

where we used a double-index notation A = (l,m) and where {•, •} is 
the Poisson bracket 

{f,g} ■= - {dtp jd^g - d^Jd^g) , p = sin<p 1 . (1) 

The factor p is conventional at this stage however it turns out to be 
a unique one (up to a constant) for which the regularization of the 
of a membrane Hamiltonian can performed. By analyzing membrane 
equations of motion one finds that, up to a constant, p = V det g rs 
where g rs is a metric of the embedded surface. 

Other reason why that factor is preferred is due to the identity for 
the coordinates of the unit sphere 

xi = cos(/?i, X2 = sin if \ cos tp2, X3 = sin tp\ sin</?2 

i.e. under a Poisson bracket defined as in ([1]) the coordinates of a 
unit sphere satisfy relations similar to those for the spin matrices. 



3 



This observation turns out to be crucial for matrix regularization of 
S 2 . Note also that the Poisson bracket satisfies all the axioms of a 
Lie bracket (i.e. the antisymmetry and the Jacobi identity) hence it 
provides the Lie algebraic structure for functions on S 2 . 

Matrix regularization maps spherical harmonics Ya to matrices 
in such a way that when replacing the Poisson bracket {•, •} by a com- 
mutator — £[■,■] the structure constants f^ c given by —i[Ta,Tb} = 

^abc^c converge in the large N limit to cabc- Let us give more 
details for that construction for S 2 . 

First step is to use the fact that spherical harmonics can be ex- 
pressed in terms of homogeneous polynomials in coordinates x\, X2, 

x 3 

^m(0>V) = ^ a ijk x l x 2 x 3 = ^a\...ai x a x ■ ■ -X ai 

i+j+k=l a;=l,2,3 

where in the last equality we introduced the coefficients a; which 
are totally symmetric. The second step is to replace the coordinates 
Xi by N x N spin matrices J, = Si/ ' \J s(s + 1), [Si, Sj] = ieijkSk where 
s is the spin (in particular 2s = N — 1 and s(s + 1) = (N 2 — l)/4 is 
an eigenvalue of the Casimir operator S 2 ) obtaining a map 

Y L^ T L'-=llN ^ C ^...a;^l • • • J CLl ( 2 ) 
Oi=l,2,3 

with the constant 7^ to be determined (see below). The choice of 
spin matrices Jj is motivated by an algebraic identity Yli Jf = 1 
which mimics the corresponding identity for the coordinates £^ • x 2 = 1 
(however one can use other conventions e.g. Jj = Si as in original work 
by Hoppe PQ). At this stage there is no reason yet to think that this 
choice will give rise to structure constants which in the large N limit 
converge to the structure constants of APD's of S 2 . Therefore at this 
stage it would be premature to refer to the above construction as to 
the construction of a matrix-regularized S 2 . 

The key observation is to note that the Poisson bracket on Xj's 
has the same structure as when considering spin matrices with the 
commutator 

£ijk x ki 

2 

- i[Ji, Jj] = -j====e ijk Jk (3) 
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(which was mentioned above) and to observe that both the Poisson 
bracket and the commutator satisfy the Leibniz rule 

{F, GH} = {F, G}H + F{G, H}, 

— i[F, GH] = —i[F, G]H — iF[G, H] (4) 

for functions F, G, H and their matrix counterparts F, G, H. Because 
of Q and (HJ), when calculating the structure constants of TVs by 
evaluating [Ta,Tb] one will perform the same type of manipulations 
as if one was calculating {Ya,Yb} modulo factors of = ~ 1/JV 

for large N. It is the main reason why the structure constants /abc 
converge to cabc- What remains to calculate is the constant 7/7V in 
((2D- It must be fixed so the all the factors of 1/N will be canceled when 
calculating an arbitrary commutation relation [Ta,Tb]- The correct 
choice turns out to be 7/7V oc N for large N with the proportionality 
factor depending on the conventions (see [23j for a detailed discussion). 



3 General considerations 

To generalize the construction outlined in the previous section, to 
the case of higher-dimensional manifolds, is a challenging task. Let 
us point out some of the most important difficulties related to such 
construction. 



3.1 The N d scaling 

First, let us note that the number of N 2 — 1 real parameters of the 
SU(N) matrix is closely related to the topological dimension of the 
surface i.e. for (i-dimensional manifold we expect the SU (N) matrices 
to be replaced by objects parametrized by ~ N d real numbers, for 
large N. The simplest way to see this is to consider a d-dimensional 
torus T d . Functions on T d can be expanded in terms of a convenient 
basis 

e w ...e' m , <^e[0,2tt) (5) 

where n^'s are integers. Because non of the (flkS are distinguished we 
expect that the regularization based on the truncation of the mode 
numbers should be such that < N for some N and all /c's. 
Therefore there would be N d parameters corresponding to a single 
expansion in terms of ([5]). 
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A similar conclusion is obtained when analyzing the expansion on 
S d in terms of corresponding S^-spherical harmonics ^ (cp\, . . . , 
These functions can be classified by homogeneous polynomials of or- 
der I of d + 1 variables [241 125] . The number of such polynomials for 
fixed I is given by 

(e.g. 21 + 1 for d = 2, (I + l) 2 for d = 3) i.e. it scales with j 3 ^y l d ~ 1 
for large I. Therefore the sum of the first N modes, for large N, gives 

N-l N-l 

Another way to argue the scaling N d is to view the regularized 
(i-manifold as a lattice of interacting points (i.e. vertices and links). 
Here one immediately arrives at the conclusion that there should be 
about N d such points. However this argument is less rigorous then 
the previous one since already for surfaces it is not clear what is a 
geometrical meaning of the N 2 — 1 parameters. Moreover as observed 
by Nicolai and Helling |22| the matrix regularization is significantly 
different from lattice approaches as it does not introduce any dimen- 
sional parameters (unlike in lattice approaches where one introduces 
a lattice spacing). 

Finally we note that the scaling N d (for large N) satisfies also 
consistency checks e.g. for S 2 x S 2 one should obtain ~ TV 4 real 
parameters, as for S 4 or T 4 , which is indeed the case - there are N 2 
modes corresponding to each S 2 in S 2 x S 2 . 



3.2 d-arrays 

The above considerations are simple but suggest something less trivial 
i.e. that the possible algebraic objects corresponding to functions on 
c?-manifolds should be certain generalizations of square matrices. A 
natural proposal is to consider d-dimensional matrices or (i-arrays i.e. 
matrices that carry more then two indices which range from 1 to N 
and a related concept of d-algebras (for a review see |26| ) . There 
exist however elementary problems with such objects e.g. for (2n + 1)- 
arrays it is difficult to define a product of two arrays. If a product 
of two arrays involves only index contractions then it is not possible 
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to multiply two (2n + l)-arrays in such a way that the result is again 
the (2n + l)-array. For example a product of two 3-arrays can be a 
4-array a 2-array or a 0-array 

(F • G)abcd ■= FabeGecd, (F • G)ab '■= FacdGdcb, 

(F»G) := FabcGcba 

but never a 3-array (here • indicates a product of d-arrays, F and 
G are d-arrays corresponding to functions F and G respectively). To 
overcome this difficulty one considers a multiplication rule in which 
some indices are fixed and equal but not summed over \27\ [28] . Similar 
problem does not appear for 2n-arrays i.e. it is possible to define a 
product of two 2n-arrays so that the product involves only index con- 
tractions and the result is again a 2n-array. For 4-arrays an example 
is 

(F • G)abcd ■= FabefGfecd 

and clearly there are many other examples of this sort obtained by 
permuting the indices and considering linear combinations. 

The apparent dichotomy between odd and even dimensions is the 
reason why regularization of d-manifolds should be carried out in 
a conceptually different manner for these two cases respectively. It 
seems also that the 2n-dimensional manifolds should be easier to deal 
with. 

Similar distinction between odd/even dimensions takes place in 
topological properties of Lie groups i.e. Betti numbers of Lie group 
manifolds are given by Betti numbers of certain products of odd 
spheres [29J. Because the structure underlying the regularization of 
manifolds is expected to be a Lie group it seems plausible that the Lie 
group cohomology is related to that distinction. 

Let us also note that in case of 2n dimensions the scaling N 2n can 
be also obtained by considering N n x N n matrices instead of In- arrays. 
There exist a natural operation, the Kronecker product of matrices, 
which gives such matrices from the usual N x N ones. Therefore there 
exists a possibility to regularize the 2n-manifolds by means of usual 
matrices which are obtained by a Kronecker product of n, N x N ones. 
We shall use this observation in Section 4. 

3.3 d-brackets 

For (i-manifolds the counterpart of membrane APD algebra is the 
cTVPD (d-volume-preserving diffeomorphism) algebra. The counter- 
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part of the Poisson bracket is the Nambu <i-bracket |30| 

{F 1 , ... , F d } := U^d Vi F x . . . d v F d (8) 

p i a 

where (fi is the manifold's parameterization, p is a certain scalar den- 
sity to be determined (in this paper we choose p as in the case of 
surfaces i.e. p is the square-root of the determinant of the metric 
of embedded manifold). Technically this is no longer an algebra in 
a usual sense (as the Nambu d-bracket defines a product among d 
functions and not just 2 functions) but a different algebraic structure 
called the d-algebra. A natural way to proceed now is to replace the 
Nambu d-bracket by an antisymmetric product of d-arrays (modulo 
some numerical factors e.g. the factor of —i in {•, •} — > ■] for the 
case of surfaces) 

{F u ...,F d }^ c d [Fx, ...,F d }:= c d e^ F h • 4 . . . • F ld (9) 

where c d is that factor. 

Let Ya be an orthonormal basis of functions (f M YaYbP = $ab) 
on a d-manifold Ai and consider an expansion of some function F on 
that manifold 

F = Y j <^aYa. 

A 

Although we do not write the range of A explicitly (due to a technical 
difficulty - the index A is a multi index and is different for different 
manifolds), the above sum is over all the modes Ya- Let us assume 
now that we have a d-array like objects T^'s corresponding to Y^'s 
hence a counterpart of F is 

F -> F:=^a A T A 
A 

where the prime indicates that the sum is over a finite number of 
modes A = A(A^) (A cx N d for large N). In order to show that such 
regularization is indeed correct one has o verify that the structure 
constants CA 1 ...A d+1 of the dVPD <i-algebra, given by 

{Ya X t ■ ■ ,YA d } = CA 1 ...A d A d+1 YA d+1 

can be approximated by the structure constants Ad+i of the reg- 
ularized (f-algebra 

Cd[T Al T Ad ) = f A ^... AdAd+1 T A d+1 



i.e that f Al !.. Ad+1 ^ c Al ... Ad+1 - 

In the following section we shall prove that assertion for particular 
4-manifolds. 

4 Regularization of 4-manifolds 

In terms of finding the appropriate regularization for 4-manifolds, the 
products of two 2-manifolds such as S 2 x S 2 or S 2 x T 2 are the simplest 
ones considering the fact that one knows how to regularize surfaces. 
From this perspective the manifold S 1 x S 3 seems to be much more 
difficult to deal with. Because the algebra behind surfaces is SU(N) 
it is reasonable to expect that the algebra underlying regularized 4- 
manifolds contains two copies of SU (N) as in e.g. the tensor product 
SU(N) (g> SU(N). On the other hand, this statement seem highly 
nontrivial for manifolds such as S 1 x S 3 . 

4.1 The Kronecker product 

Let us start with S 2 x S 2 . We shall embed this manifold in R 6 and 
parametrize its coordinates (x±, X2, X3, yx, 7J2, 2/3) as 

xi = cos(/?i, X2 = sin^i cosv9 2 , X3 = sin tpi sin tp2 , 

yi = cos ip 3 , y 2 = sin^ 3 cosv7 4 , y 3 = sin ip 3 sin c/? 4 . (10) 

Functions on this manifold can be expanded in terms of products of 
spherical harmonics for each sphere 

^(^1,^2,^3, ¥>4) = ^2 a MA 2 YA 1 { l Pl-, t P2)YA 2 {^?„ l PA) 
Ai,A 2 

where we again use the double index notation A\ = (Zi,mi), A2 = 
(h, m 2)- If we now replace Y A 's with matrices T A as in ([2]) then a 
question arises what matrix operation is the counterpart of the prod- 
uct Y Al ((pi,(p2)Y A2 ((p3,(p4). A usual matrix product i.e. 

Y Al {vi,tp2)Y A2 (ip 3 ,<p A ) -)• T Al -T M 

is not good. The reason is simple - T Al T A2 is still N x N matrix 
hence the number of independent parameters describing T Al T A2 scales 
with iV 2 while we require the iV 4 scaling. However the result of the 
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Kronecker product of two matrices (A B)^^ := AijB^i has the 
correct scaling, therefore we consider 

^ (^1,^2)}% (923,^4) -» lhi 2 NT Al ®T A2 (11) 

where 7y 2 jv is a constant to be determined later on (a counterpart of 
1lN in ©)• That this prescription may work follows form bi- linearity 
of the Kronecker product and the mixed-product property 

(A x (8) A 2 ) ■ (Bi <g> B 2 ) = (Ax ■ Bx) (g> (A 2 ■ B 2 ) (12) 

which, together with the fact that TVs are given by linear combination 
of products of spin matrices Jj, implies that we can expand T Al ®T Ai 
in terms of (ordinary) products of matrices 

Jf := Ji ® 1, Jf := 1 <g) Ji- 

Concretely, we have 



T Al ®T M = ( £ C. a;i J ai . . . J a;i ® £ c™ 2 ..^ J 6l . . . J l 

,^=1,2,3 / \6;=1,2,3 



^2 C ™i-a h J ai • • • J a h I ' ^..^ J bi " " " J h 2 ) ■ ( 13 ) 

^=1,2,3 J \bj=l,2,3 

4.2 The double-Poisson bracket 

Note that the matrices Jf and Jf form two independent sets of N 2 x 
./V 2 spin matrices i.e. 

r jx jxi _ 2i jx rjy T yi _ 2i T y 

[Jf,Jj]=0 (14) 

which implies that there is no ordering ambiguity in (|13p . Moreover let 
us observe that the commutator •] corresponds now to a double- 
Poisson bracket i.e. defining 

if, 9} = {/,5W 2 +{/,5W4! {f,g}a/3 = -^-(d a fd p g - dpfd a g) 

(15) 
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we obtain relations on S 2 x S 2 which are similar to (|14p 

{xi,Vj} = 0. (16) 

Therefore any further construction must take into account the fact 
that the commutator ■] is a regularization of the double-Poisson 
bracket (fT5|) . 

4.3 4-arrays 

The advantage of writing everything as in (|13p is that we have got rid 
of the explicit Kronecker product by " moving" it inside the definition 
of the building blocks Jf , Jf. These matrices can now be manipulated 
in a usual fashion however since they are iV 2 x N 2 the resulting reg- 
ularization will be much different from that for surfaces. This result 
is a realization of what was discussed in section 3.2: one can view the 
Kronecker product / ® g as a 4-array with the indices ranging form 1 
to N 

or as a 2-array (a matrix) with indices from 1 to N 2 

(f®g)u, i,J = i,...,n 2 

(we use the notation F or f®g for matrices corresponding to functions 

F(Vl,<P2,<P3,<P4,) or f(tpx,tp 2 )g(ip3,ip4,) on S 2 x S 2 ). 

4.4 General prescription 

Up to now we have defined a map 

Xi ->■ Ji ® 1, yi ->• 1 ® Ji (17) 

which gives us the matrix counterparts of the basis Ya 1 Ya 2 via (fTTI) 
and (|13p . However for Cartesian product of arbitrary 2-manifolds this 
prescription will not work. On the other hand, we could have obtained 
the same result by simply defining Ya 1 — > T J 4 1 ® 1, Ya 2 — > 1 <8> Ta 2 
and by requiring that the product of functions YajYa 2 is mapped to 
a (usual) product of iV 2 x N 2 matrices (T Al ® 1)(1 ® Ta 2 ). This 
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prescription will work since (T Al <8> 1)(1®T A2 ) = T J 4 1 8) T Aa and hence 
can be used in the general case. 

Let us therefore consider two 2-manifolds N and M and their 
product N x M.. Functions on M and M. can be expanded in terms of 
modes Y^j (</?i, (/J2) and Y A2 {ip 3: ^4) respectively. Suppose their matrix 
counterparts are T J 4 1 and Ta 2 and define a map 

Y Al (<pi,(p2)->T Al ®l, Y A2 (ip 3 ,ip 4 )^l®f A2 . (18) 

Having done that we now define the counterpart of the Nambu 4- 
bracket by 

{^,^2,^3,^4} -> c 4 [F 1 ,F 2 ,F 3 ,F 4 }:=c 4 e» kl F i F j F k F l (19) 

where i^'s are the matrix counterparts of i^'s and C4 is a number to be 
determined. The product FiFjFkFi in (fl~9|) is a usual product of four 
TV 2 x N 2 matrices. In the following we will be interested in calculating 
the Nambu 4-bracket for the basis functions Y Al Y A2 and their matrix 
counterparts T Al (g) T A2 therefore all Fi 's considered here will have that 
product form (and all F^s will have the Kronecker product form). 

We now use the observation that the Nambu 4-bracket can be 
expanded as [3l] 

{Fi,F 2 ,F 3 , F4} 

= {F 1: F 2 }{F 3 , F^ - {F 1 ,F 3 }{F 2 ,F 4 } + {F u F 4 }{F 2 , F 3 } (20) 

where {•,•} is a double-Poisson bracket, and that the quantum 4- 
bracket (the 4-commutator) resolves as 

[Fi, F2, F 3 , F4] 

= {[A, F 2 ] , [F 3 , F 4 ] } - { [A , #3] , [h ih]} + { [F l , F 4 ] , [F 2 ,F 3 ]} (21) 

i.e. it has twice as many terms due to the noncommutativity. There- 
fore it is reasonable to expect that —\ [•,-,-,•] is a correct quantum 
counterpart of the Nambu 4-bracket (with the minus factor due to 
the % factor in {•, •} —> —*[•,•] and the fact that the 4-commutator is 
a sum of "squares" of commutators, implying that C4 = —1/2, c.p. 
(fl9l) ). For the case of S 2 x S 2 we have shown that the double-Poisson 
bracket corresponds to the commutator — £[•, ■]. Let us verify that this 
rule is true for general product manifolds. Considering the functions 
in the product form and their matrix counterparts 

F\ : = /lO/^^siO^,^) -> fi ®gi =■ A, 
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we find that the double-Poisson bracket and the commutator are 
{F\,F 2 } = {fx, f2}ipnp 2 9l92 + f2fi{gi,g2}<p 3 <p 4 , 

-i[F 1} F 2 ] = (-i[hJ2})®gi92 + hfi®(-i[9i,92}) 

i.e. they resolve in the same way (note the ordering of fi's). This, 
together with (|20|) and (|2ip . verifies that the Nambu 4-bracket can be 
replaced by — 5 [ v > v] while keeping in mind that the double-Poisson 
bracket is replaced by •]■ 

4.5 Leibniz rules 

Let us now note that the Leibniz rule for a double-Poisson bracket i.e. 

{F 1 ,F 2 F 3 } = {fi, f 2} Vltp2 f 3919293 + M/1,/3}^ 

+ hhh{gi, 92} mVi 93 + hhh92{gi:93]ip^A- (22) 
resolves in the same way as the commutator 

-i[F 1 ,F 2 F 3 ] = {-i[fi, / 2 ]/3) ® hfafa + h{-i[fi, M) ® 3iff2#3 

+ /2/3/1 ® (-*[&,&])& + /2/3/1 ® &(-*[&,&]) (23) 

(where F 3 = / 3 (pi, ^2)53(^3, ^4) A ® ff3 := ^3). Eqs. ((22]) and 
((23|) can be used in evaluating {-^1,-^21-^3,-^4-^5} and its quantum 
counterpart — F 2 , -F3, -F4-F5] via identities ((20]) and (f2T]) respec- 
tively. As a result one obtains a Leibniz-like rule, analogous to dU) 
for surfaces, which can be used to evaluate the Nambu 4-bracket 
and the 4-commutator for basis functions Ya 1 Ya 2 and their matrix 
counterparts Ta 1 ® Ta 2 respectively. Just like in the case of S 2 we 
observe that when performing that calculation one makes the same 
algebraic manipulations modulo factors of 1/N (for large N) com- 
ing from the commutators [fi, fj] and [gi,gj] (c.p. (fT4l) for the case 
of S 2 x S 2 ). Therefore in the leading order in TV the structure con- 
stants f^iAiAsAiAs approximate the 4VPD structure constants 
CAiA 2 A 3 A 4 A 5 - The large N expression for the factor 7; 1 z 2 iv is given by 
the product Jifj^jf^ i-e. of order iV 2 . 

Note that identities ([22]) and ([23]) are useful here only because the 
basis functions on N x M and their matrix counterpart are in the prod- 
uct form and Kronecker product form respectively. However the basis 
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corresponding to e.g. S 4 will not have that property and therefore 
one has to use/argue a different mechanism which contains Leibniz- 
like rule. Unfortunately the Leibniz rule for a Nambu 4-bracket i.e. 

{Fx, F 2 , F 3 , -F4F5} = {Fi, F2, F 3 , F^F^ + F 4 {Fx,F 2 , F 3 ,F^} 

does not hold for the 4-commutator. On the other hand, from matrix- 
regularization point of view one may relax the Leibniz rule to allow 
1/N corrections - i.e. the following identity would still work 

[Fi,F 2 , F 3 , = [A,F 2 , F 3 ,F 4 ]F 5 + F 4 [Fx,F 2 ,F 3 ,F 5 ] + 6 (24) 

where O are some -Fj-dependent terms which vanish in the large N 
limit. Using (|2ip and the Leibniz rule for the commutators we find 
that matrix O is 

6 = [Fx,F 4 ][F 5 , [F 2 ,F 3 \] + [F 3 ,F 4 ][F 5 , [F lt F 2 ]] - [F 2 ,F 4 ][F 5 , [F U F 3 ]] 

+[[F 1 ,F 2 \,Ft][F 3 ,F s ] - [[Fi,F 3 ],F 4 ][F 2 ,F 5 ] + [[F 2 ,F 3 },F 4 }[Fx,F 5 ]. 

All the terms in O contain three commutators hence they are of order 
1/N 3 (assuming that each commutator is of order 1/iV) while the 
4-commutators is of order 1/N 2 hence we verify that O is indeed 
subleading in N. 

Therefore in calculating the structure constants one can use the 
Leibniz rule for 4-commutator keeping in mind that there will be ad- 
ditional 1/N terms. This is another way of proving that the structure 
constants of 4VPD are recovered when N — > 00. 

4.6 S A 

To be more specific and show that the finite N Leibniz rule (|24l) is 
required, let us consider a unit S embedded in R 5 and given by co- 
ordinates 

xx = cos tpx, x 2 = sin (fx cos ip 2 , x 3 = sin^i sin<^2 cos ip 3 , 
X4 = sin tpx sin ip 2 sin ip 3 cos (p 4 , x$ = sin tpx sin tp 2 sin ip 3 sin ip 4 . 

Functions on S can be expanded in terms of S^-hyperspherical har- 
monics y^ im2m3 (tpx, ¥2, tp 3 , tfi). These harmonics can be expressed in 
terms of homogeneous symmetric polynomials in variables Xx, x 2 , x 3 , X4 
and X5 as 

vl — n mim,2m s i 1 i 2 is 14^15 

1 m 1 m 2 m 3 ~ u ii 1213*415 1 2 ^3 x 4 x 5 

U+«2+i3+«4+l5=^ 
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mim 2 m 3 



a fc = l, 2,3,4,5 

where the coefficients c™ 1 ™^" 13 are symmetric in indices a^'s. Because 
the coordinates xi for S do not separate as in (fTU|) a map similar to 
(|17p will not work. On the other hand, the Nambu 4-bracket can be 
evaluated for S 4 coordinates, we have 

{%i ; Xj , Xfc , X[\ — tijklm-Em j ^, J, A?, / , 771 — 1,2,3,4,5 

where the scalar density used here is p = sin 3 ip\ sin 2 if 2 sin ^3 (c.p. 
©)• Therefore a matrix counterpart of Xj, call it Tj, should satisfy 

~~ 2 ^j; Tfc) T/] = d&€ijki m T m (25) 

where g?4 may depend on N (dj should scale like 1/N 2 based on ar- 
guments from previous subsection). Here we have assumed that the 
formula for the quantum Nambu 4-bracket should be independent of 
the topology of the manifold (as it is in the case of surfaces) therefore 
we use (|19p with C4 = —1/2. An obvious choice of matrices satisfy- 
ing (j25j) is Ti oc 7j where 7j's are 5D Dirac matrices corresponding 
to Euclidean signature {7i,7j} = 2<5yl - identity ([25]) follows directly 
from algebraic properties of 74 's. However Dirac matrices in 5D are 
4x4 and there are no other irreducible representations of a different 
dimension while we are looking for N 2 x iV 2 matrices Tj for any N. 
It follows that matrices Tj cannot satisfy Clifford algebra for N > 2 
while for N = 2 they are given by 5D gamma matrices. 

The last remark suggests the following construction, define 

r * = a (Jk f)' r4 = 6 (i 0)' r5 = 6 (o -1 

(26) 

where Si's are n x n spin matrices [Si,Sj] = ze^^Sfc, a and 6 are 
normalization factors so that = 1- These matrices are 2n x 2n 

hence to obtain the correct N 4 scaling one should take 2n oc N 2 for 
large N. To argue again the importance of the iV 4 scaling let us also 
note that since there will be about 2N /4! matrix counterparts of the 
harmonics l£ lfnama (c.p. ©) 

Y l —>T l —'Sim r m i m 2 m 3r r ('271 

- 1 raim2m3 r - L mim2rri3 ■ IU\ / ^a\...ai a-i • • • a-i V ' I 

o fe =l, 2,3,4,5 
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the majority of them would be linearly dependent if e.g. 2n ~ N. The 
scaling 2n ~ N 2 guaranties that they are independent. 

For 2n = 2, these matrices are proportional to the 5D Dirac ma- 
trices in the chiral representation, however for 2n > 2 they no longer 
satisfy Clifford algebra (e.g. {ri,r 2 } 7^ for 2n > 2). On the other 
hand, as we shall show now, Tj's satisfy identity (|25p with 

I 6 1 6 / 6 

a -yW^Ty 71' d '-^5^Ty (28) 

For the proof one should a priori calculate twenty four 4-commutators 
in (|25p however due to the complete antisymmetry it is enough to 
evaluate only 5 of them. Using identity (|21h we find that they are 

[Ti, r 2 , r 3 , r 4 ] = -4a 3 s ( s + i)r 5 , [r 1; r 2 , r 3 , r 5 ] = 4a 3 s(s + i)r 4 , 

[ri,r 2 ,r 4 ,r 5 ] = -i2a6 2 r 3 , [ri,r 3 ,r 4 ,r 5 ] = i2a6 2 r 2 , 
[r 2 ,r 3 ,r 4 ,r 5 ] = -i2a6 2 r! 

where s(s + 1) is an eigenvalue of n x n matrix representation of the 
Casimir operator S 2 and 2s = n — 1. The normalization condition 
^T, r| = 1 and the fact that d 4 in (125p should be common to all the 
4-commutators imply equations 

a 2 s(s + 1) + 26 2 = 1, -12a6 2 = -ia 3 s(s + 1) = -2d 4 

which are solved by (|28p (note that other branches are possible e.g. 
a — > —a, b — > —b, d 4 — > — d 4 ). Therefore we have found a 2n x 2n 
irreducible matrix representation of the 4-algebra (|25p . As anticipated 
the factor cf 4 scales like 1/n ~ \/N 2 therefore in the large N limit the 
4-commutator is zero which is analogous to the case of S 2 (c.p. ([3])) i.e. 
as N — > 00 matrices Tj 4-commute. They will not commute however 
since 

i r - r »i=!(S "o 1 ) < 29 > 

which holds for all n. 

Having defined matrix counterparts of Y^ im2m3 by (f2"7j) we can now 
evaluate corresponding structure constants. Once again the scaling 
2n ~ iV 2 turns out to be necessary - otherwise most of the structure 
constants would be zero (for 2n ~ N) or there would be too many 
degrees of freedom (for 2n ~ N k , k > 2). One can now use the 
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Leibniz rule with 1/N corrections (|24[j 1 l to prove that the structure 
constants given by T4 im2m3 's converge to the structure constants of 
the 4VPD of S 4 . Finally let us observe that form ([6]) it follows that 
there are 

N-l 

^2D l4 = -N(N + l) 2 (N + 2) 
1=0 

matrices XJri im2TO3 which for large N is less then (N 2 — l) 2 hence the 
group underlying S 4 is not SU(N) <g> SU(N). 

5 S 3 

Although in this paper we are mostly interested in 4-manifolds let 
us observe that matrices (|26p can be used to regularize S 3 . For the 
coordinates of a unit 5 3 embedded in M 4 we take 

x\ = cosy?!, X2 = sin 991 cos 9?2 5 £3 = s\mp\ sint/?2 cos (^3, 
X4 = sin ifi sin </?2 sin </?3 

which satisfy 

{xi,Xj,x k } = -eijkixu k, I = 1,2,3,4 (30) 

with the scalar density p = sin 2 (pi sin ip2 ■ Just like in the case of S 4 , 
in order to find a matrix counterpart of f|30|) it is clear to start with 
the 4D Dirac matrices j{ owing to the identity [7i,7j,7fe] oc tijul<oli- 
However because of the factor of 75 one cannot identify Xj's with 7j's 
yet. A way to proceed further is to alter the definition of the quantum 
3-bracket and define it using the 4-commutator as 

[F,G,H] 15 :=[F,G,H, l5 ] (31) 

which therefore gives [7^, 7j, 7fc] 75 oc tijuili and hence Xi can be re- 
placed by 7i's. This approach was used in the context of multiple- 
membrane theories [32] where it is also shown that it arises naturally 

1 Since the commutator |29|) does not scale like 1/N one may think that the matrix O in 
([Ml) for a generic case may be of order 1/N 2 (instead of 1/N 3 ) and hence not subleading. 
However O will in fact scale like 1/n 2 ~ 1/N 4 due to r a 's, a — 1,2, 3, carrying the factor 
of 1/n. 
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from non- associative algebras. Let us now generalize this construction 
to 2n x 2n matrices. Using the same IYs as for S we find that 

[Fi) r.j) rfc]r 5 := [Ti, Tj, Fk, F5] = — Id^ijkSmFm = ZdiEijkiFi 

therefore one may use these matrices for negative (i.e. the branch 
a — > —a, b — > —b, — > —d^) to represent quantum counterparts of 
coordinates of S 3 as long as we are using [•, •, -]r 5 as a definition of a 
quantum 3-bracket. We observe that in the large n limit not only the 
[Fi, Fj, rfc]r 5 's but also the commutators [1^,1^] are 0. 

The spherical harmonics of S 3 can now be mapped to 2n x 2n 
matrices via. 

Y L im2 (vu ¥2, <pa) = c ^.™( x «i • • • 

o fe =l, 2,3,4 

^T l mim2 -.= llN Y, C. m 4r ai ...r a; (32) 

Ofc=l, 2,3,4, 

and therefore to maintain the iV 3 scaling of the number of modes we 
will take 2n ~ N 2 . As for S 14 , matrices Fi are A^ 2 x iV 2 but the 
resulting matrix-harmonics <\32h contain only A 3 degrees of freedom 
(more precisely, it follows form ([6|) that there are A(A + l)(2A + l)/6 
degrees of freedom corresponding to first A modes). 

That the structure constants coming from T4 im , 2 's converge to 
the structure constants of VPD of S* 3 follows immediately from the 
fact that the quantum 3-bracket is in fact given by a 4-commutator. 
Therefore the reasoning presented in the case of S 4 applies here. 

6 Summary and Outlook 

Matrix regularization of embedded surfaces is a procedure that at first 
sight seems easily generalizable to ti-dimensional manifolds. One sim- 
ply performs Fourier expansion of functions on a manifold and then 
replaces the modes by a suitably chosen matrices - why should it be 
more complicated? The problem however is there - it is much more 
difficult to find concrete representations of algebraic structures (dif- 
ferent from Lie algebras) that should be used to replace the Nambu 
(i-bracket. Even more basic problem is related to noncommutative ob- 
jects that replace functions on the manifold. These objects should de- 
pend on roughly N d real parameters where A is the number of Fourier 
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modes - a natural choice would be to use the d-arrays. Although we 
find this path appealing, in this paper we focused on possible ma- 
trix representations due to the ambiguity in defining the product of 
two d-arrays (d > 2). For 2n-manifolds we simply consider N n x N n 
matrices while for d = 2n — 1 a subset of N n x N n matrices. 

In this paper we focused on d = 4 manifolds due to the possible 
relevance of this case to quantum gravity. By finding a matrix regu- 
larization of space-time one may in principle be able to quantize the 
resulting system of finite degrees of freedom. We firstly considered 
manifolds that are Cartesian product of two 2-manifolds and showed 
that functions on such manifolds can be approximated by N 2 x N 2 
matrices from SU(N) (g) SU(N). This result shows that one does not 
have to refer to 4-arrays in order to regularize 4-manifolds. We then 
considered more complicated case: S" 4 , which also can be regularized 
by iV 2 x iV 2 matrices however this time we were unable to identify 
the underlying group. Finally we used the construction for S 4 to find 
the regularized description of S 3 . In particular we used the observa- 
tion that one can define a quantum 3-algebra by using the quantum 
4-bracket (the 4-commutator). 

Let us now discuss some generalizations. First, as for the product 
manifolds M X M for arbitrary closed manifolds A/*, M of dimension 
> 1 we conjecture that the resulting group of matrices underlying the 
regularization is Gj^f® G\m> where the groups Gj\f, Gm correspond to 
M and Ai respectively. For the proof, the construction presented in 
Section 4.4 will work provided the corresponding Leibniz rule (with 
1/N corrections) is also there. 

Second, the regularization of S , d > 4 seems straightforward by 
appropriate generalization of matrices T, (|26p . The strategy is to first 
regularize S with the usual 2n-commutator as a quantum counter- 
part of the Nambu 2n-bracket and then apply this construction for 
S 2n ~ l . For example S 6 embedded in M. 7 can be regularized by taking 
An x 4n matrices Tj (i = 1, . . . , 7) with the Nambu 6-bracket replaced 
by the 6-commutator. After this is done the regularization of S 5 can 
be made by using the same matrices (excluding Tj) and defining the 
quantum 5-bracket by [•,-,-,•, -]p 7 = [•,-,-,•, •,!*•?]. 

Third, as observed in [28] there exist 3-array, n x nxn representa- 
tions of the 3-bracket for S 3 . That construction has a virtue of using 
the usual definition of the quantum 3-bracket (i.e. the 3-commutator). 
Therefore one does not need to refer to the 4-commutator [-,-,-,75]. 
On the other hand, working with 3-arrays requires introducing a non- 
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standard multiplication rule. It seems plausible that these two ap- 
proaches are equivalent. 
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